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Abstract 

We propose a hybrid model of simplicial quantum gravity by performing at once dynamical 
triangulations and Regge calculus. A motive for the hybridization is to give a dynamical descrip- 
tion of topology-changing processes of Euclidean spacetime. In addition, lattice diffeomorphisms 
as invariance of the simplicial geometry are generated by certain elementary moves in the model. 
We attempt also a lattice-theoretic derivation of the black hole entropy using the symmetry. Fur- 
thermore, numerical simulations of 3D pure gravity are carried out, exhibiting a large hysteresis 
between two phases. We also measure geometric properties of Euclidean 'time slice' based on a 
geodesic distance, resulting in a fractal structure in the strong-coupling phase. Our hybrid model 
not only reproduces numerical results consistent with those of dynamical triangulations and Regge 
calculus, but also opens a possibility of studying quantum black hole physics on the lattice. 
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1 Introduction 



One of the most difficult problems in modern physics is the construction of a consistent theory of 
quantum gravity, whereas Einstein's theory of classical gravity has been very successful in explaining 
the large scale structure of spacetime PP. Many programs to formulate the full quantum theory of 
gravity are under active research Since none of them are decisive at present, it is still an open 

problem to understand what is quantum spacetime at microscopic scales. 

In this paper we explore a possibility that dynamical Regge calculus 01^1, which is a hybrid lattice 
model of dynamical triangulations [3 IBJ and quantum Regge calculus [HI , gives a possible candidate for 
a constructive definition of quantum gravity. Although traditional approaches, namely, Regge calculus 
and dynamical triangulations, have been well studied for a long time as lattice field theories of gravity, 
they are not satisfactory in several respects. A necessary enlargement of physical degrees of freedom 
in lattice gravity strongly motivates our study on the hybridization. Actually, the enlargement enables 
us to define exact "diffeomorphism-invariance" on the lattice at least classically 

We will apply the symmetry to a lattice-theoretic derivation [S] of the Bekenstein-Hawking entropy 
of a black hole. By the end of 1970's it was generally accepted that a black hole has the huge entropy 
5bh proportional to its horizon area A [17\ \1H\ : 

'^^^ = I^^t^ ' 

where is Boltzmann's constant, c the speed of light and G Newton's constant^. Nowadays, it is 
generally believed that if one sticks to usual theories of gravity (e.g. Einstein's gravity and dilaton type 
gravity coupled to normal matter fields), eq. is widely valid 22|. Theoretical explanation of the 
origin of the huge entropy is considered to be a necessary ingredient for any consistent theory of 
quantum gravity pUj . 

Another motive for studying the hybrid model is giving a description of topology-changing processes 
of Euclidean spacetime in a dynamical way. In a continuum approach Hawking calculated semiclas- 
sically the contribution of gravitational instantons to such processes of four-dimensional Euclidean 
spacetime manifolds and discussed the phenomenon in connection with Regge calculus; the resul- 
tant spacetime that is highly curved and has all possible topologies is called the spacetime foam. Such 
a complicated spacetime is expected to appear in the strong-coupling region of quantum gravity 
especially inside black holes and in the very early universe. Thus, a question naturally arises: How can 
we describe the spacetime foam in lattice quantum gravity? It will be shown that the topology-changing 
processes can in principle occur via degenerate simplicial configurations in our hybrid model. 

Regge's lattice formulation of general relativity |^ has been applied to quantum gravity in two 

different manners, that is, quantum Regge calculus and dynamical triangulations, as mentioned above. 

In the former, all link- lengths in a fixed pattern of triangulation (a fixed connectivity of vertices) 

play the role of dynamical variables instead of the metric field g^^. Accordingly, the integration of 

the link-lengths with a proper functional measure is assumed to give a constructive definition of the 

^We write down explicitly all the constants in eq. Q in order to show how large the black hole entropy is. Indeed, for 
a black hole of horizon area A — 1 cm^, one obtains a huge value Ssn/kB ~ lO''^. 
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path integral for the field g^i,. In the latter, in contrast, all patterns of triangulations are regarded as 
dynamical, while all the link-lengths are fixed to a single lattice spacing a. In this case, the sum over 
all possible triangulations is expected to give another constructive definition of the path integral for 
the gravitational field. 

Both lattice models of gravity have merits and demerits; the weakest point common to both of them 
is the lack of the gauge symmetry that corresponds to 'reparametrization-invariance' on the lattice. 
Actually, one must first give a precise definition of the 'reparametrization-invariance' on the lattice in 
order to construct a 'gauge-invariant' measure. We call the symmetry lattice dijfeomorphism-invariance. 
Physically, the invariance should be the lattice counterpart of the principle of general covariance. An 
intention of this paper is to discuss the lattice symmetry even on the finite lattice by covering as large 
degrees of freedom as possible in the space of Riemannian geometries. Thus, one would conjecture that 
such an enlargement of degrees should be realized by performing at once the link-integration arising 
from quantum Regge calculus and the triangulation-sum arising from dynamical triangulations. 

The plan of this paper is as follows. In section [2 we give a short review of quantum Regge calculus 
and dynamical triangulations, making a brief comparison between them. How each formulation regu- 
larizes the space of Riemannian geometries is the central issue there. Section 01 deals with the definition 
of the hybrid model, that is, dynamical Regge calculus, and discuss its fundamental properties. It will 
also be discussed that the topology change of discrete spacetime can occur in principle via degenerate 
simplicial complexes, because the lattice action for the gravitational field remains well-defined and finite 
on any degenerate configuration that is not a manifold. Section 0] considers the symmetric property, 
that is, the lattice diffeomorphism-invariance. We will also attempt a lattice-theoretic derivation of 
the black hole entropy using the lattice symmetry and give a simple interpretation of it in con- 
nection with the spacetime uncertainty principle of string theory |49[ I50| I51j . Section El is devoted 
to numerical studies of three-dimensional pure gravity. In particular, a fractal structure based on a 
geodesic distance (31] is measured in detail. As a result, we will obtain a picture of quantum spacetime, 
which is similar to that obtained in DT. Furthermore, we will acquire several other numerical results 
consistent with those of both dynamical triangulations and quantum Regge calculus. In section El we 
will summarize our results and discuss several perspectives of dynamical Regge calculus. 



2 Quantum Regge calculus vs. dynamical triangulations 

In this section we will make a comparison between the two formulations of lattice gravity. We first 
explain briefiy the idea of the Euclidean path integral approach in the continuum, and then present a 
basics of Regge calculus. Emphasis is laid on geometric properties on the finite lattice. 

2.1 Regge's formulation 

Leaving aside the lattice approaches for a moment, let us look at the Euclidean path-integral 
approach in the continuum The basic object of the approach is the functional integral for the 
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metric field g^i, of Euclidean signature + + •••+: 

where Seh^] is the Einstein-Hilbert action on a d-dimensional Euclidean spacetime M: 

Here yl is the cosmological constant, and units are such that c = h = 1. For simplicity, we assume M 
to be a compact, closed manifold of a fixed topology. The integration ((2)) is taken over the space of the 
metric field g^i, with an appropriate functional measure "Dg. 

The groundwork for the path integral approach was laid by De Witt [17j and, subsequently, de- 
tailed prescriptions were developed by Gibbons, Hawking and Perry ^|. Although this approach has 
several difficulties and unsolved problems, a number of highly suggestive results have been obtained. 
Presumably, the dramatic successes of this approach are that, as mentioned in introduction, the foam 
picture of spacetime can be established semiclassically by calculating the contribution of the gravita- 
tional instantons 53 , and that the creation of particles near a Schwarzschild black hole can be related 
in a direct and simple manner to the properties of the Euclidean Schwarzschild solution |19j . 




Figure 1: An example of the simplicial decomposition. A smooth manifold M with a metric g is 
replaced by a simplicial manifold T with a set of link-lengths {li,l2, ■ ■ ■ ,k, ■ ■ ■ , l-Ni}- is the dehcit 
angle around a hinge (vertex) h, corresponding to the scalar curvature. In two dimensions the volume, 
Ah, of the hinge h is dehned to be 1, while A^ takes non-trivial values in higher dimensions. 

Having shortly surveyed the continuum approach, let us now return to our main interest. In 

Regge's lattice formulation of general relativity 30_, the continuum spacetime M is replaced by a 

discretized space T that consists of a finite number of d-simplices, as shown in Fig. ^ Such a discrete 

space T is called a d-dimensional simplicial manifold or triangulation in lattice gravity^. First, we 

give a connectivity of vertices in T, which must satisfy the manifold condition that the space looks 

locally like a Euclidean space R'^. Incidentally, the total number, Nk, of fc-simplices also is determined 

(k = 0, 1, . . . , d). Next, we give lengths li, I2, ■ ■ ■ , Ini to all the links (1-simplices) in T. As a result, 

^Exactly speaking, the discrete counterpart of a Riemannian manifold {M,g^,j) is a piecewise linear (PL) manifold, 
which has a certain metric structure that converges on the Riemannian structure in a continuum hmit |16| . In other 
words, a simphcial manifold equipped with the PL metric structure is assumed to be a (discrete) physical spacetime. 
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the pair of the connectivity and the set of all the link- lengths determine the simplicial geometry 
on the discrete spacetime T. This procedure is called the simplicial decomposition. An example of the 
procedure is depicted on Fig. where a two-dimensional torus with a metric g^u is decomposed to a 
simplicial torus with the link-lengths {/j} by gluing triangles (2-simplices) along their boundary links 
(1-simplices). 

On the simplicial manifold, in general, one can define geometric quantities (differential forms, 
curvatures, parallel transport, etc.) in a standard way. In particular, the Einstein-Hilbert action Q is 
replaced with the so-called Regge action, 5'Regge > on T [HOj : 



where is the volume of a {d — 2)-simplex (hinge) h and the volume of a d-simplex a. The 
dimensionless quantity €h in the right hand side of (jlj) is called the deficit angle around the hinge 
h, which plays the role of the local scalar curvature on T. The symbol I denotes the set of all the 
link-lengths {/i, h, ■ ■ ■ , Ini}- 



Figure 2: An example of a simplicial complex C that is not a manifold. The complex C is a union 
of three simplicial submanifolds Ti [i = 1,2,3), and obviously violates the manifold condition at the 
junction parts, namely, at the vertex a and along the link be. The Regge action Sneggei^, Ti] is well- 
defined and finite on each Ti. Furthermore, the Regge action S'Reggef^ C] on C is also defined uniquely 
as the sum of the actions J2i 'S'Reggc[^ Ti]. 

Interestingly, the Regge action Q is well-defined even on a degenerate configuration that is not a 
simplicial manifold but a simplicial complex. The crucial difference between the simplicial manifold 
and the simplicial complex is that the former satisfies the manifold condition that the space looks 
locally like a Euclidean space R^, but the latter does not necessarily. Fig. [21 shows a two-dimensional 
simplicial complex C which is not a manifold of definite dimension; the manifold condition is violated 
at the junction parts. However, the Regge action S'Reggel^ C] on C is uniquely defined as the sum of the 
actions 'S'Reggel^, Ti] where Ti {i = 1, 2, 3) are the simplicial submanifolds of the complex C. Unlike 
the case of the continuum action Q, the finiteness of the Regge action (@J even on such a degenerate 
complex is a notable feature, particularly for describing the topology change of Euclidean spacetime in 
later section 




(4) 
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2.2 Quantum Regge calculus 

From now on let us concentrate on the problem of lattice quantization of gravity based on Regge 
calculus (RC). In quantum RC the discretized equivalent of the integration over the metrics Q is 
implemented by varying the link-lengths on a fixed triangulation , assuming that patterns of triangula- 
tions are not dynamical at all; this is the reason why quantum RC is often called the fixed triangulation 
(FT) approach. Accordingly, the partition function for quantum RC on a fixed triangulation T is given 



where dfixi^ is a link- integration measure defined on T, including constraints imposed by the (higher- 
dimensional analogs of) triangle inequalities. Each length li is integrated over a range /min < k < ^max- 
In eq. © we explicitly write down the 'subscript' T of the measure dfij^m and the 'variable' T for the 
partition function Zrc[^] in order to make it clear that they are defined on the triangulation T. This 
notation is useful for our formulation. 

The FT approach © has been actively applied to quantization of gravity ^ ^] and actually 
has many appealing aspects. However, this approach has some subtle issues. One of them is how 
to define the discretized measure dfj,T[l] that corresponds to the continuum measure "Dg in eq. (|2|) 
over the so-called superspace of the Riemannian metrics on M. In fact, one has to gauge away 
the diffeomorphism group DiffM that is the gauge freedom of general relativity. In the continuum a 
possibility is to start from the supermetric formulation [3^ that defines a gauge-invariant measure over 
the superspace and then to gauge away DiSm by using one's favorite gauge. 

Here one would ask a simple question: What is the lattice counterpart of the diffeomorphism 
group Diff^^? Though we have no clear answer to the question, the notion of "diffeomorphism" on 
the lattice has been used at two different levels ^37. One definition of the symmetry, which is often 
adopted by those studying Regge calculus in classical relativity, is called "invariance of the geometry" 
in which transformations of link-lengths leave the geometry invariant. A proper implementation of 
the definition is to require that all local curvatures (and hence all deficit angles) should be unchanged 
under the transformations. But this requirement is too strong to satisfy in the FT approach^ |12j . 

Another definition, favored by those wishing to use results from lattice gauge theories, is referred 
to as "invariance of the action" in which transformations of link-lengths leave the action invariant. It 
is possible even in quantum RC to imagine that changes in the link-lengths which could increase deficit 
angles in one region and compensatingly decrease them in another would produce no overall change in 
the Regge action (jlj. However, the symmetry in this sense holds approximately at most in (almost) 
flat space ^DEIE]; still less it would be possible to flnd in curved space any exact symmetry that 
can be interpreted as the lattice diffeomorphism in this approach. 

The other problem is that there can be generally singular conflgurations including such simplices 

that have very short and very long links at the same time, as shown in Fig. |H| Though one usually 

■^The only exception is flat space where an infinite number of choices of the link-lengths will generally correspond to 
the same flat geometry even in the FT approach. 
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deals with this problem by putting upper and lower bounds on the lengths li, such configurations 
are unavoidable to obtain a metric with very high curvature out of the flat metric. It is still unclear 
whether one can explore regions in the space of metrics where the metric is very singular and fairly 
different from the typical metric on the reference simplicial space that we have chosen at the beginning. 
It will cause a large difficulty in studying the strong-coupling region of lattice quantum gravity, where 
configurations with very high curvatures are expected to emerge frequently. 




Figure 3: Singular configurations arising in the FT approach, (a) A singular 2-simplex. (b) A singular 
3-simplex. 

2.3 Dynamical triangulations 

In the dynamical triangulation (DT) approach [7| |5j, which is the alternative to the FT one, it 
is assumed that none of the link-lengths in each triangulation are dynamical and all of them can be 
fixed to a single lattice spacing a. Accordingly, the simplicial manifold T considered in DT consists 
of equilateral d-simplices. The Regge action (jlj) on such a equilateral triangulation T becomes the 
following simple form: 

SBT[T] = -KoNo + fidNd , (6) 

where kq and are constants related to G, A and a. Nq and stand for the number of 0-simplices 
and that of d-simplices on T. Instead of the link-lengths, we take varying connectivity of vertices 
as dynamical. Hence, the path integral for DT is defined by the sum over all possible patterns of 
triangulations that consist of equilateral d-simplices: 

Zdt = ^ exp(^-Sr,T[T]^ ■ (7) 

Titriangulations 

Here we fix the topology of all the triangulations T. 

In what follows, we make a brief comparison between FT Q and DT ((T)) to clarify the differences 
between them. For any Riemannian manifold non-singular enough to start as the reference at the 
beginning, one can construct a sequence of simplicial manifolds. We expect that in DT the discrete set 
of simplicial spaces is regularly distributed in the space of all the Riemannian metrics |33j , as shown in 
Fig. 13(a). In contrast, simplicial spaces arising in FT are localized in the neighborhood of the smooth 
metric on the regular lattice that has been chosen at the beginning (see Fig. |l](b)); such a localization 
will make the lattice configurations less accessible to the strong-coupling regions that might be far from 
the reference metric, and prevent "diffeomorphism-invariance" from holding. In this respect, DT has 
the advantage over FT. Moreover, DT has the natural UV cutoff a of the theory. 
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(a) (b) 




Figure 4: Sketch of the space of aU Riemannian metrics. Each point means a simpUcial geometry. 
The circles with dot represent the reference (Hat) simpUcial space at the beginning, (a) "Regular" 
distribution arising in the DT approach, (b) "Localized " distribution around the reference metric 
arising in the FT approach. 

Nevertheless, the DT approach makes some issues more difficult. Actually, we have completely lost 
diffeomorphism-invariance on the lattice, because one cannot deform smoothly (even approximately) 
a simplicial lattice into another in this approach. Some argue that in DT the permutation group 
acting on Nq vertices is expected to reproduce the diffeomorphism group acting on a smooth 
manifold in a continuum limit^, although no exact proof has been known yet. One has to check 
whether diffeomorphism-invariance is recovered in the continuum limit^, if it exists. In addition, it is 
not clear how to define the classical continuum limit of DT, since this approach gives no field equation 
corresponding to Einstein's equation. As a consequence, this difficulty makes it fairly intractable to 
study (quantum) black hole physics in DT. 

Given reasonable geometric and symmetric properties, universality in quantum field theory will 
generally assure the same results for physical observables in a certain continuum limit, even though 
the details of UV behaviors of two field theories are different from each other. This is almost the case 
with lattice field theories, where physical results are expected to be independent of the specific details 
of the UV cutoff and the specific forms of the lattice actions as long as theories of interest preserve 
the underlying symmetries. In the case of lattice gravity, however, such a naive expectation from 
universality will not necessarily hold, because the two approaches have different symmetric properties 
as discussed above. Exactly speaking, it is a drawback to the lattice regularizations of gravity that 
they lack diffeomorphism-invariance even at the classical level, in contrast to other lattice field theories 
where gauge symmetries are (classically) exact on the finite lattice. Attempts to overcome the drawback 
strongly motivate us to study the hybrid model of quantum RC and DT, namely, dynamical Regge 
calculus, from now on. 

*In the IIB matrix model, which is a candidate for a constructive definition of superstring theory, the same scenario is 
also discussed |42|. 

^In the case of two-dimensional quantum gravity, however, it has been verified that results obtained in DT ^ coincide 
with the predictions from CFT |25| and diffeomorphism-invariance is recovered in two dimensions. 
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3 Dynamical Regge calculus 



In this section we first give the definition of our hybrid model as a simultaneous implementation 
of quantum RC and DT and discuss its fundamental properties. We discuss possible behaviors of 
the entropy® of the model, which is related directly to the well-definedness of the model. Then, we 
formulate elementary local moves necessary to carry out numerical simulations. Finally, we will give an 
extension of the hybrid model including degenerate configurations and apply it to a topology-changing 
process of Euclidean spacetime on the lattice. 

3.1 Definition of the hybrid model 

The partition function, Zdrc, of dynamical Regge calculus (DRC) is defined by performing at once 
the link- integration © and the triangulation-sum Q: 

Z-DRC = X] / ^I^tII] exp(^-S'Regge[^, T^j = ^ Z^iclT] , (8) 

T:triangulations r:triangulations 

where each link-length Li is integrated over the range Imin < < ^max- The sum X/'j':trianguiations 
means that any pattern of triangulation (connectivity of vertices) T with various sets of the link- 
lengths is included so long as it satisfies both the manifold condition and the triangle inequalities (see 
Fig. [51) • The reason for the name of dynamical Regge calculus is that Regge calculus Zrc [T] on each 
triangulation T is thought of as if a 'dynamical variable' in the partition sum (jHI- 




Figure 5: Lattice configurations appearing in the sum Q. The Erst two configurations have differ- 
ent patterns of connectivity, arising from the triangulation-sum ^rp. The next two have the same 
connectivity but contain different lengths of the i-th link, arising from the link-integration J dfiTil-]- 

A few remarks are in order. First, the behavior of the entropy of the hybrid model (jSJ is essential 

to its well-definedness as a statistical system. A natural generalization from DT leads to the following 

®Here the entropy means the total number of possible triangulations in the model and, therefore, it is not the black 
hole entropy. Do not be confused with the same term 'entropy'. 
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definition of the entropy, W{Nd), of the model 

T:Na fixed 

where the sum ^T-Na&xcd taken over triangulations with the number of d-simpHces A'',^ fixed. It 
depends strongly on the behavior of the entropy Q whether the model is well-defined statistical- 
mechanically or not. If we set the coupling constant G = oo for simplicity, the following inequality 
holds: 

where the minimal length Zmin plays the role of the UV cutoff. Thus, for the partition function ZDRci^^d) 
with N[i fixed, one easily obtains the following inequality: 

T:Nd fixed 

< W{Nd) e-^'-in^'' . (11) 

If an exponential bound for the entropy W^Nd) holds 

W{Nd) < const. X e'^'^'min^d (yi^ . ^ positive constant) , (12) 

then DRC (jH)) is well-defined at least as a statistical system. In this case it is a possibility that one 
can take a continuum limit by fine-tuning the cosmological constant A. However, it is very hard to 
give any analytic proof for the exponential bound H12j) because the integration measure dfiT[l] includes 
the triangle inequalities which are too intractable to calculate analytically. Instead, we will later see a 
piece of numerical evidence that the bound (fT2|) holds in case that we use the scale-invariant measure 

Second, we comment on a direct relation between DT and DRC ((SJ. If one chooses the 5-function 
measure 

di-iT[l] = Yldli6{li-a) , (13) 

i 

then DRC Q becomes the same as DT (O where all the lengths are fixed to the single lattice spacing 
a and the triangle inequalities are automatically satisfied. In other words, DT is equivalent to DRC 
with the (5- function measure (|13() . As is well known with matrix models j26j, such a choice (|13() gives 
a constructive definition of two-dimensional quantum gravity, even though the measure (|13jl breaks 
explicitly the lattice diffeomorphism-invariance. Although the problem of the lattice symmetry is not 
dealt with here, we will observe in section |2 that it is a difficult task to construct a 'diffeomorphism- 
invariant' measure on the lattice. 

3.2 Hybrid (p, q) moves in DRC 

Here we construct local, ergodic moves, i.e., local changes of triangulations, though they never 
change the topology of the simplicial manifolds; such elementary moves are a necessary ingredient to 
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carry out numerical studies of DRC (jH)). The Monte-Carlo method is applicable to DRC"^, 'S', as well 
as it has been to both quantum RC 113j and DT ^Sj- Remember that in DT the so-called {p, q) 
moves^ are used and they are ergodic in the class of triangulations of fixed topology [301 ■ what 
follows, we give an extension of the {p, q) moves to invoke quantum RC in addition to DT according 
to DRC ©. 




arbitrary lengths 
(b) vo I 




Figure 6: Hybrid {p, q) moves in two dimensions, (a) A hybrid (1, 3) move, where the triangle voViV2 
is divided into three triangles vvqVi, vviV2 and vv2Vo. Then, inserted links vvq, vvi and VV2 can take 
arbitrary lengths so long as triangle inequalities hold. Its inverse move, called the hybrid (3, 1) move, 
is also depicted, (b) A hybrid (2, 2) move acting on two triangles vqViv^ -\- V1V2V3. After deleting the 
link V1V3 of length I, a new link vqV2 of length I is inserted. The length I is also arbitrary unless triangle 
inequalities are broken. 

In Fig. ini we show an example of the extended moves in two-dimensions; in Fig. IHl (a), a triangle 

VQV1V2 is divided to three new triangles vvqVi, vviV2 and vv2V() by inserting a vertex v and three new 

links vvq, vvi, vv2. The link-lengths are arbitrary as long as triangle inequalities hold. The deficit 

angles around the vertices vq, vi, V2 and v can change continuously in DRC as well as in quantum 

RC^. The inverse move is defined straightforwardly, as shown in Fig. [HI (a). We call the two moves 

of Fig. El (a) the hybrid (1, 3) move and the hybrid (3, 1) move. Similarly, one can extend the (2, 2) 

move of DT to a hybrid one of DRC, as depicted on Fig. |H1 (b). The link ^1^3 of length / shared by two 

triangles ^0^11^3 and V1V2V3 is flipped to a new link V0V2 of length /; / also is arbitrary unless triangle 

^In ref. j40|, these moves are called the (fc, I) moves. To avoid confusion, however, we use here the term of the {p, q) 
moves instead, since we use the symbol I to denote (a set of) link-lengths. 

^Indeed, we observe soon later that the link-integration of quantum RC is automatically included in such moves. 



11 



inequalities are broken and, hence, the deficit angles around the vertices vo, vi, V2 and can take 
continuous values. We call the extended move in Fig. (b) the hybrid (2, 2) move. 

It is straightforward to generalize the two-dimensional hybrid {p, q) moves to higher-dimensional 
ones, where p + g = d + 2 (d = 3, 4). For example. Fig. [7| shows three-dimensional hybrid (p, q) moves. 




arbitrary lengths 



I 




Figure 7: Hybrid {p, q) moves in three dimensions, (a) A hybrid (1,4) move, where a new vertex v is 
put into the 3-simpIex vv()ViV2V3 and four hnks vvq, vvi, vv2 and vv^ with arbitrary lengths Iq, li, I2 
and I3 are inserted unless triangle inequalities are broken. Its inverse move is also depicted, (b) A 
hybrid (2, 3) move acting on two 3-simplices vqViV2V^ + viV2V^V4. Instead of the triangle viV2V^, the 
new link vqv^ is inserted. Length I of vqV4 is arbitrary so long as triangle inequalities hold. Its inverse 
moves, called the hybrid (3, 2) move, is also shown. 

The hybrid {p, q) moves described above are expected to be ergodic in the sense that all the lattice 
configurations can be generated only by the hybrid {p, q) moves. Moreover, in Fig. |H1 we explain how 
the hybrid moves implement not only DT but also quantum RC . A link-integration will be done by 
successive applications of the hybrid moves as follows: 

Ml) On a two-dimensional simplicial lattice we pick up a link uv of length li (see FiglH] (a)). Then, 
we apply a hybrid (2, 2) move to the link uv, and a new link ab is created (see FiglHl (b)). 

M2) Another hybrid (2, 2) move is applied to the link ab, getting back to the same connectivity as the 
initial configuration (see Fig|Hl(c)). The only difference between the initial and final configurations 
is that the length li of the link uv is replaced with the different one I / . 
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Through the processes Ml) and M2), the length of the hnk uv has changed from the initial li to the final 
one If according to the proper Boltzmann weight. Hence, this is a typical example of the realization 
of quantum RC only by a combinations of the hybrid moves. 

u 




Figure 8: An example of the realization of a link-integration by successive applications of the hybrid 
{p, q) moves in two dimensions. Through the first hybrid (2, 2) move Ml), the link uv of length li is 
hipped to a link ah. The next hybrid (2, 2) move M2) changes the link ah back to the link uv of length 
If. As a result, this process realizes a link-integration of quantum RC only by the hybrid moves. 

One would guess that in higher dimensions similar combinations of the hybrid moves can incorporate 
link-integrations of quantum RC into DRC. This is really the case with our hybrid model. In addition 
to the link- integration discussed above, the hybrid {p, q) moves are ergodic in the sense that all the 
patterns of triangulations can be generated by the hybrid moves owing to the established ergodicity of 
the usual {p, q) moves in DT jlU]. As a consequence, these properties complete the ergodicity of the 
hybrid {p, q) moves in DRC. 

3.3 Description of topology change in DRC 

As an application of the hybrid model, we try to describe the topology-changing processes of Eu- 
clidean spacetime on the lattice. For the purpose, we need to extend our model © in order to deal 
with degenerate configurations that appear in the topology change, as will be discussed below. 

Over forty years ago. Wheeler pointed out [HH that the Einstein-Hilbert action (j^J allows large 
fluctuations of the metric and even of the topology of spacetime manifolds on scales of order of the 
Planck length. This is due to the fact that the action for the gravitational field (jSJ is not scale 
invariant, unlike that for the Yang-Mills fields. Hence, a large fluctuation of the metric over a short 
length scale does not have a very large value of the action © and so is not highly damped in the 
path integral The resultant quantum spacetime with the large fluctuations of both the metric and 
the topology is called the spacetime foam. Subsequently, Hawking calculated semiclassically the path 
integral for the spacetime foam by summing up "gravitational instantons^" with various Euler numbers 
X in four dimensions 53 . Physically, the "gravitational instantons" describe the topology-change of 
Euclidean spacetime. Furthermore, Hawking discussed a dynamical realization of such a topology- 
changing process by using quantum Regge calculus; a metric can change topology without increasing 
the lattice action by more than an arbitrary small amount [HS]. In what follows, we explain the close 
relation between Hawking's idea and our hybrid model. 

^For a detailed explanation of the gravitational instantons, see ref. |54|. 
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We first define the partition function, .^edrc> of the extended hybrid model: 

^EDRC = X] / '^'^'^t^l exp(^-5Regge[/, C]^ = ^ ^RCp] , (14) 
Cicomplexes Cxomplexes 

where the sum X^c compiexes / dfJ-cil] admits degenerate simplicial complexes (see Fig. EJ in addition to 
simplicial manifolds. Furthermore, manifolds (or complexes) of various topologies are contained in the 
extended sum (|14)) in order to describe topology-changing processes dynamically. 

Now we discuss a simple example in Fig.Elto illustrate a process of changing continuously from the 
topology of one simplicial manifold to the topology of another. Consider a simplicial manifold S'^ with 
a 'waist' triangle abc (Fig. (a)). The three links ab,ac and be of the 'waist' has lengths labJac and 
Ibc respectively, satisfying a triangle inequality lab < lac + he- The manifold has Euler number x = 2, 
where x is defined on the two-dimensional lattice as 

X = No-Ni + N2 . 

Here stands for the number of /c-simplices {k = 0,1,2). Once an equality lab = lac + he holds 
by varying the link-lengths, the triangle abc will collapse to a singular link along which the manifold 
condition does not hold any longer (Fig.|Ul(b)). In general, if some of the (sub)simplices collapse to lower 
dimensions, a simplicial complex will not remain a manifold but become a degenerate configuration; 
this is the case with the simplicial complex shown in Fig. 121 (b). 

Subsequently, we delete four 2-simplices adc, bdc, bee and eac around c, while combining the two 
links ae and eb into a new link ab. Then, c is deleted. Moreover, a separation of the new link ab 
(dashed) from the old ab (solid) will lead to other degenerate complex shown in Fig. 121 (c) . The shaded 
part is an empty region on which no (sub)simplex resides. Although the simplicial complex is no longer 
a manifold, the Regge action @ remains well-defined and finite even on such degenerate configurations 
shown in Figs.|Hl(b) and (c). The process from (b) to (c) is beyond the link-integrations (of RC) and 
the hybrid {p, q) moves (of DRC); in general, such a process should be called a surgery acting on the 
simplicial complex. 

Next step is to apply hybrid (p, q) moves several times to the simplices around a and b. As a result, 
we can obtain a simplicial complex as shown in Fig. |^ (d), in which coordination numbers of a and b 
are three. However, the obtained complex still has apparently Euler number x = 2, though it is not a 
manifold of definite dimension and topology. 

Finally, we "blow up" the degenerate simplicial complex to obtain a new simplicial manifold. We 
delete six 2-simplices around a, that is, aaia2, 00203, 00301, aa[a2, oOgOg and oogo'^ in Fig. (d). 
Then a is deleted, while we identify Oj with a[ (i = 1, 2, 3) pairwise. Similarly, we delete six 2-simplices 
around 6, and identify bi with 6^ (i = 1,2,3). As a consequence of the "blowing-up" surgery, the 
degenerate complex of Fig. (d) has changed to a simplicial manifold of Fig. |H1 (e) . 

In this way, one can pass continuously from one metric topology to another with the Regge action 
remaining finite. Indeed, one could deform the topology of the simplicial manifold only by the action 
of such local moves within the framework of extended DRC (|14|) . However, we meet with a difficulty 
in finding a minimal set of local moves that are enough to generate step by step all the (degenerate) 
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(a) manifold 5^ (x = 2) (b) complex {x = 2) 




Figure 9: Example of a topology-changing process via degenerate complexes in extended DRC p4l) . 
(a) A simplicial manifold S'^ of x = -^as a 'waist' triangle abc. (b) If three links of abc satisfy 
lab = lac + he, then abc collapses to a (singular) link. The configuration is now not a manifold but a 
degenerate complex, (c) The four 2-simplices sharing the 0-simplex c are deleted. But we leave the 
two 1-simplices ac and be and combine them into another 1-simplex ab (solid curve), and then delete c. 
The dashed curve is the old ab. The shaded region between the two ab becomes empty, (d) By acting 
several hybrid (p, q) moves around a and b, the two vertices will have coordination number three. The 
complex still has x = 2. (e) We delete six 2-simplices around a and further six 2-simplices around b. 
Then, the 0-simplices ai and a[ (i = 1, 2, 3) are identified and, similarly, bi and b[ (i = 1, 2, 3) are done. 
The complex has changed to a simplicial manifold T^ of x = by the "blowing-up" surgery. 
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complexes of various topologies. This is due to the surgery operations, such as the process from 
Fig- El (b) to (c) or the process from Fig.|Hl(d) to (e). Therefore, it is a challenging task for us to study 
numerically the extended model (|14j) . 



4 Lattice diffeomorphisms in DRC 

Given a Riemannian manifold (M, g), there is arbitrariness in choosing coordinate-systems that 
cover the manifold M. The arbitrariness is necessary to represent the principle of general covariance 
in general relativity. Similarly, in the case of lattice gravity, the freedom would be reflected as the exis- 
tence of the infinite number of triangulations that correspond to the same Riemannian geometry; they 
should be transformed to each other by "lattice diffeomorphisms" defined on the simplicial manifold. 
An important issue in regularized theories of quantum gravity is the nature of the symmetric properties. 
Though significant discussions on the issue have been given, no systematic formulations of the symme- 
try have been established. Actually, in the FT approach diffeomorphism-invariance is approximately 
realized as the appearance of gauge zero modes only on (almost) flat backgrounds ^lElE]; by the 
method of lattice weak-field expansion, the zero modes and the corresponding eigenvectors identified 
with infinitesimal local coordinate transformations in a continuum limit. But invariance properties on 
curved backgrounds are beyond the scope of the weak-field expansion. 

In this section we will see that the Regge action ^ on arbitrary curved backgrounds is exactly 
invariant under certain hybrid {p, q) moves of (extended) DRC and further that such moves are inter- 
preted as the lattice diffeomorphisms in a natural way. 

4.1 Invariance {p, q) moves on the simplicial complex 

As an illustrative example, we first construct a hybrid (2, 2) move under which the Regge action Q 
remains exactly invariant. We zoom up some 2-simplices around a hinge (vertex) /i on a triangulated 
surface, as shown in Fig. IIUI We put them on a flat plane to make explicitly visible the deficit angle 
efi around h, as depicted on Fig.^l(a). A hybrid (2, 2) move makes the link hj flip to other link ik in 
such a way that the vertices i and k are connected by a straight line on R^, as shown in Fig. (b). 
Evidently, the deficit angles eh, e^, ej, et around h, i, j, k are invariant, because the following relations 
hold: 

ef + ef = 9f , = + ef) . (15) 

where oi"^^ and oij"^ {v = h,i,j,k; m = 1,2) are shown in Fig. (a) and (b). The sum of areas of 
the triangles also is invariant: 

Vhij + Vhkj = Vhik + Vijk ■ (16) 

Here the volumes (a = hij, hkj, hik, ijk) denote areas of the 2-simplices considered. These simple 
relations (|15() and H16|) ensure that the local geometry and the Regge action @ are exactly invariant 
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Figure 10: Example of the invariance (2, 2) move which keeps the simphcial geometry and the Regge 
action (0) invariant, (a) Some triangles sharing the hinge h are put on a Bat plane in order to 
visualize the dehcit angle eh- (b) The link hj is hipped to other link ik in such a way that the length 
of the link ik is the Euclidean distance between the vertices i and k. 

under the move of Fig. IIUI Hence, we call this hybrid (2, 2) move an invariance (2, 2) move. From 
a geometric viewpoint, the invariance (2, 2) move is a simplicial representation of a local coordinate 
transformation, because the 2-simplices hij, hjk, ihk and ijk themselves can be regarded as local 
coordinates and the move changes one local coordinate system spanned by hij and hjk to another 
spanned by ihk and ijk. Furthermore, it keeps exactly invariant the simplicial geometry in the simple 
way. Therefore, we can naturally interpret the invariance (2, 2) move as a lattice diffeomorphism in 
the sense of "invariance of the geometry" . 




Figure 11: Example of the invariance (1,3) move, which keeps invariant the simplicial geometry and 
the Regge action 0). (a) A triangle ijk is put on a fiat plane, (b) The triangle is subdivided into 
three triangles hij, hjk and hki in such a way that the length of the link hi is the Euclidean distance 
(straight line) between h and i and, in the similar way, the lengths of hj and hk are chosen to be the 
Euclidean distances. 

Similarly, we can define other invariance moves as shown in Fig. 1111 We first pick up a 2-simplex 
ijk shown in Fig. ^2(a), and then subdivide it into three simplices hij, hjk and hki in Fig. ^2(b) by 
putting a new vertex h inside ijk and connecting it with i, j, k in such a way that the flatness property 
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inside ijk is preserved. Under the action of this move, the invariance of the deficit angles e^, ej, 
around i, j, k is guaranteed simply by the following relations: 

0, = ^f) + 0f), ej = 0f^ + 0f\ ^fc = ^f + 0f. (17) 

Here, the dihedral angles 9y and oij^'' {v = i,j,k; m = 1,2,3) are shown in Fig. (a) and (b). The 
deficit angle eh around the inserted vertex h is exactly zero, because the Euclidean flat geometry inside 
the triangle ijk is preserved exactly, namely, the following equality holds: 

27r = ^i') + + . (18) 

The total volume also is invariant under the move: 

Vijk = yhij + Vhjk + yhki , (19) 

where the volumes denote areas of the 2-simplices a = ijk, hij, hjk and hki in Fig. 1111 We call 
this hybrid (1, 3) move an invariance (1, 3) move, and its inverse an invariance (3, 1) move. These 
relations (|17|) . H18() and ()19() guarantees the exact invariance of the local geometry under the moves. 
These moves can be interpreted as the lattice diffeomorphisms in the same way as the invariance (2, 2) 
move could be. 

It is straightforward to generalize the two-dimensional invariance {p, q) moves to higher-dimensional 
ones, under which both the simplicial geometry and the Regge action are kept exactly invariant in the 
similar way. For example, we consider the three-dimensional case. Let us remember Fig. where the 
3-simplex V0V1V2V3 is divided into the four 3-simplices by the hybrid (1, 4) move. If one adjusts the 
lengths Iq, li, I2 and ^3 so that the geometry inside VQV1V2V2, remains flat, the deficit angles around the 
six links of vqViV2V^ are invariant and, furthermore, deficit angles around the four inserted links are 
also zero (flat) in Fig. [7| (b). The sum of volumes of the four new 3-simplices is exactly equal to that 
of vqViV2V^. This action is an invariance (1, 4) move that keeps the Regge action @ invariant in three 
dimensions. Similarly, we can construct an invariance (2, 3) and an invariance (3, 2) moves. 

In addition to the symmetry described above, the permutation group, S^^, of Nq vertices gives 
another invariance, corresponding to the re-labeling of vertices on the lattice. Therefore, we conclude 
that the lattice diffeomorphisms in the sense of invariance of the geometry are generated by both the 
invariance {p, q) moves and the permutation group S^q in DRC. Symbolically, we denote the symmetry 
of DRC as follows: 

00 

|lattice diffeomorphisms I = |moves generated by invariance {p, q) moves| k S^o ■ (20) 

What does this symmetry imply to lattice quantum gravity? One might naively expect that one can 
simply integrate over the gauge transformations (the lattice diffeomorphisms) without any gauge-fixing 
procedure just as one could in lattice gauge theories. However, this is not the case with gravitation |llj : 
the crucial difference between Einstein's gravity and gauge theories is that the symmetry group is 
nan- compact in theories of (lattice) gravity. Thereby, one must factor out the infinite volume of 
the diffeomorphism group to make sense of the path integral for the gravitational field. Thus, it is 
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plausible that the infinite volume of the lattice diffeomorphisms would prevent the entropy bound (|12() 
from holding. Generally speaking, it is very difficult to prove (or deny) the inequality ()12|) by analytic 
methods, because the entropy © depends on the triangle inequalities which are too complicated to 
estimate analytically. Moreover, it will also depend on the measure chosen for the link-integration. 
Hitherto, we deform the problem to another one: what measure enables the exponential bound 1)12(1 to 
hold? In section 121 we will investigate numerically the (un)boundedness of the entropy © using a few 
types of measures. 

4.2 Black hole entropy in DRC 

Before turning to numerical studies, let us now apply the lattice symmetry to the black hole entropy 
problem. Among various attempts to solve the problem, Carlip advocated that by virtue of the horizon 
structure, the "would-be gauge freedom" of general covariance supplies physical states of the Kerr 
black holes in arbitrary dimensions [SSI- spite of some shortcomings in his analysis [201 ? idea 
has attracted much attention to the relation between the entropy and the horizontal geometry of black 
holes. Indeed, it is argued in ref. j57j that the smooth diffeomorphism on the event horizon can be 
regarded as a nontrivial asymptotic isometry of the Schwarzschild black hole. In what follows, we 
discuss a simple method of deriving the black hole entropy from a lattice-theoretic viewpoint; the 
derivation is based on the lattice diffeomorphism ()2U() in DRC. 

Suppose that there exists a simplicial lattice that corresponds to a four-dimensional black hole with 
a triangulated horizon 5^ of area A. We define Ghm as a set of all the moves generated by the hybrid 
{p, q) moves acting on the triangulated black hole. We define a subset Hum C Ghm as 

i^HM = {-^ ^ G'hm ; / ^ {lattice diffeomorphisms}, /I52 7^ identity, and 

/I52 S {lattice diffeomorphisms on the triangulated horizon 5^ }| . (21) 

An example of an element / € Hum is depicted on Fig. ^1 Physically, each element / G Hum means 
such a quantum fluctuation that is a nontrivial combination of the hybrid {p, q) moves acting on the 
triangulated horizon; its restriction, /I52, on the horizon is required to be a two-dimensional lattice 
diffeomorphism to keep invariant the horizontal geometry. 

Assuming that there exist np elements of Hum per the Planck area ip = Gh/c^, we can easily count 
the total number of such fluctuations, Nfi^^, around the horizon in a combinatorial way: 

Nfluc ~ (np)^/'p . 

This number Nji^c corresponds to the total number of quantum fluctuations that keep invariant the 
simplicial geometry on the horizon. Hence, we can define the entropy, Sun-, of the black hole as 

Sun = ku log Nfiuc ~ ku log np ■ . (22) 

ip 

Our idea of the derivation of the black hole entropy (|22() is based on the simple picture that any 
quantum states are associated with quantum fluctuations; the classical picture of the horizon should 
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/1 52 is a lattice diffeomorphism on the horizon 




triangulated horizon S 




links on the horizon 
links around the horizon 



Figure 12: Example of a combination of the hybrid moves f £ Hhm defined by eq. i|2I|) . On the 
triangulated horizon S"^ , f\g2 is required to be a two-dimensional lattice diffeomorphism. 



accompany many fluctuations around it at a quantum level, and the huge entropy of the black hole 
arises from the contribution of such fluctuations. 

Here, we encounter a difficult problem: If one takes a continuum limit, the fluctuation density np 
will diverge, resulting in an infinite value of the entropy 5bh in eq. (|22() . Such a divergent behavior 
of the black hole entropy often appears also in continuum approaches. In particular, Susskind and 
Uglum showed that the entropy per unit area of a free scalar field propagating in a fixed black 
hole background is quadratically divergent near the horizon and such quantum corrections to the black 
hole entropy are equivalent to quantum corrections to the gravitational coupling constant G, although 
the theory of interest is non-renormalizable. They concluded that the question on the finiteness of 
the entropy is inextricably intertwined with the renormalization of the coupling G and, therefore, the 
finiteness property cannot be understood without the complete knowledge of the ultra-violet behavior 
of the theory. 

Taking the discussions in the continuum into consideration, we are led to the following idea: If 
we regard the lattice model (jH)) as an effective theory with a finite cutoff the fiuctuation density 
np will remain finite in eq. (|22|) . This finiteness enables us to avoid the divergence behavior and, 
furthermore, to normalize the Planck length ip in eq. (|22jl so as to set 

np = e^/^ (= 1.284...) . (23) 

As a consequence, we reach the following expression: 

Sbh ~ . (24) 

Thus, eq. (|24j) agrees with the Bekenstein-Hawking area law of the black hole entropy 

How can we interpret the relations ((221) and (|24() on physical ground? Interestingly, the notion 
of the space-time uncertainty principle has been advocated in the development of string theory and 
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noncommutative field theory |3^. According to the principle, space-time has an uncertain property in 
itself; uncertainty of temporal coordinate AT and that of spatial one AX satisfy 

AT-AX>£l , (25) 

where ig is the so-called string scale of order 0{ip). The lengths AT and AX, called the extremal 
length, are defined in a conformally invariant way [IHl- An interpretation of the relation (|25)) is given 
in the context of noncommutative geometry ^j, in which the coordinates are promoted to operators 
that satisfy commutation relations similar to those of quantum mechanics. However, we take another 
interpretation of the uncertainty principle (|25jl: space-time has the intrinsic minimal length^'' £s and 
we identify it with the UV cutoff /min of DRC: 

ATj-j^iYi — ^^min — ~ /min • (^6) 

As a result, the uncertainty principle can be satisfied without introducing noncommutativity of space- 
time. According to this interpretation, eq. (|23|) means that on the minimal area l'^^^ ~ £p we observe 
about one quantum fluctuation of the (lattice) gravitational field. 



5 Numerical study of 3D pure gravity 

In this section we turn to numerical calculations in order to study non-perturbatively the ground 
state of DRC (jSJ and the behavior of the entropy ©. We will calculate some observables (integrated 
curvature, average link-length, and surface area distribution function characterizing a fractal structure) 
to study three-dimensional pure gravity based on DRC. 

First, we rewrite the measure dfiTl^ in the following form convenient for our numerical study: 

/^l /-'max , s 

df^T[l] = n / dk exp(-5[/, T]j 6t{A) , (27) 
i=\ ^min 

where StI^) stands for the constraints by the (higher-dimensional analogs of) triangle inequalities, and 
we denote the measure-induced part of the action as e~'^[''^l. Two examples of the measure-induced 
part S[l, T] will soon appear. From now on we set the minimum lattice spacing /min to be one^^. 

In the DT approach the number of d-simplices, A^^, fluctuates largely in higher dimensions {d = 
3, 4), making statistics of numerical data fairly worse. Similarly, such fluctuations occur also in our 
hybrid model. In order to obtain better statistics, we add an extra term (A^ — V)^ to the Regge action 
Q); it controls the volume fluctuations around the designed size V. Hence, the total action, S[l, T], on 
the lattice is of the form: 

S[l, T] = SRegge [/, T] + S[l, T] + ^{Nd - Vf . (28) 

^"Actually, it is discussed in ref. |51| that a particular uncertainty relation leads to the minimal length for each component 
of the space-time coordinates, namely, eq. (I26II may be satisfied also in string theory. 

^^If dimensional arguments are needed, the appropriate power of the lattice spacing (which has the dimension of length) 
can always be restored at the end of calculations by invoking dimensional arguments. 
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Here 7 is a small constant. As a result, the partition function of DRC becomes 



^DRC = Yl / n^^^ exp(-S[/, T]) <5t(A) 



(29) 



T:triangulations i 



where the sum 



Titriangulations 



admits only simplicial manifolds satisfying the manifold condition. 



Numerically, the sum of lattice configurations appearing in eq. H29|) is carried out through the hybrid 
(p, q) moves, as explained in section EiH (see Fig. [7)). 

From now on we focus on three-dimensional pure gravity. The topology of the simplicial lattice is 
fixed to three-sphere and we never consider the topology-changing process here; this is due to the 
algorithmic difficulties in simulating manifolds of various topologies, which remains to be solved. 

5.1 Phase structure of 3D pure gravity 

In general, local operators cannot be gauge-invariant observables in quantum gravity. But some 
global quantities are gauge-invariant, for example, the average scalar curvature per volume: 



Here is the deficit angle around the i-th link (hinge) /j and being the volume of a 3-simplex a. 
This quantity (R) plays the role of the order parameter. The average length of all the links is also a 
good observable: 



These observables ()3U() and ()31|) can be calculated by the standard Monte-Carlo technique. In what 
follows, we perform numerical simulations using two types of measures, that is, the uniform measure 
d^j and the scale- invariant measure Y\^dli/li. The number of 3-simplices A''3 is almost fixed at 
around V = 5000, and each length /j is integrated over an interval 1 < < 10. The coefficient of the 
extra term is set at a value 7 = 1.0 x 10^^. 

5.1.1 Calculation using the uniform measure 

First, we use the uniform measure in simulating DRC (jEJ: 

dfiT[l]=lldli5TiA) . (32) 

i=l 

In this case, the lattice action S[l, T] is of the following simple form: 




(30) 




(31) 



S[l,T] = Snegge[l,T]+^{N3-Vf . 



(33) 



The partition function becomes 




(34) 
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A = 40 




5560 5570 5580 5590 5600 
number of 3-simplices 

Figure 13: The N-^ distribution using the uniform measure dli. The couphng constant G is set to be 
oo. Different curves correspond to different values of the cosmological constant ranging from A = 40 
to A = 100. The maximum value of the horizontal axis is the upper cutoff of N3. 

However, the measure ()32() causes a severe problem that the hybrid system (|34() will not be well-defined 
statistical-mechanically. This is due to the reason that the entropy bound ()12() for the configuration 
entropy Q will not hold in this case. Now we must discuss this point in more detail. 

Actually, we tried carrying out numerical calculations using the measure ()32() . In our numerical data 
with G = 00 (strong coupling limit), however, the number of 3-simplices N3 exhibits a clear tendency 
to diverge, even though we set large values of the cosmological constant A, as shown in Fig. 1131 Though 
the value of the cosmological constant varies from 40 to 100 in Fig. I13( the N3 distribution is squeezed 
to the upper cutoff = 5600. Even if we choose larger values of the cutoff, the divergence behavior 
still appears in the same way. Hence, we cannot control the system (|34j) at all. From a viewpoint of 
statistical mechanics, one can interpret such a pathological behavior as follows. If the entropy bound 
(|12() does not hold, the configuration entropy WlN^) (given by eq. @) will dominate the system in the 
strong-coupling phase (high-temperature phase). In other words, lattice configurations with large 
will frequently appear owing to the large values of the entropy Vl^(A''3), leading to the divergence of the 
lattice size. Actually, this is the case with our calculation using the uniform measure in Fig. 1131 The 
large volume of the lattice diffeomorphisms will make large (presumably infinite) the entropy of such 
pathological configurations. In that sense the data of Fig. indicates a piece of numerical evidence 
for the existence of the lattice symmetry. 

At any rate, our calculation using the uniform measure cannot proceed further. Instead of the 
pathological measure, we should next try to use other measure that may break the diffeomorphism 
invariance at the quantum level. 
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5.1.2 Calculation using the scale- invariant measure 

Next, we simulate the system H29() using the scale- invariant measure: 



dli 



i=l 



(35) 



Obviously, this measure is invariant under the local rescalings li — > Cj/j where c, are constants. In this 
case, the partition function is defined as 

"10 ^1 



yscale—inv 
^DRC 



Y,j nd/iexp(-5[/, r])5r(A) 

T ^ i=l 

where the action S[l, T] is of the form: 

S[l, T] = 5Regge[/, T]+Y, log k + l{N^ - Vf . 



(36) 



(37) 



First of all, we measure the distribution, as shown in Fig. 1141 One observes in the figure that Nj, 
distributes smoothly around the central value V = 5000, in contrast to the uniform-measure case shown 
in Fig. 1131 According to the numerical data of Fig. 1141 one can naively expect that the exponential 
bound (|12|) will hold under this measure. 
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Figure 14: A^3 distribution using the scale-invariant measure Yl^dli/li. 

In this case, we can obtain in principle any observables numerically, unlike the case of the uniform 
measure. Indeed, we measured the scalar curvature per volume, {R), as shown in Fig. 1151 Evidently, 
this system has two phases, namely, the strong and the weak coupling phases; the solid line represents 
the cooling process from the strong coupling phase to the weak one, and the dotted its inverse. One 
can clearly observe a large hysteresis in Fig. 1151 suggesting the first-order nature of the transition. It 
is consistent with results obtained in numerical studies of three-dimensional DT j43j . 
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Figure 15: Hysteresis for the curvature {R) under the scale-invariant measure Y\i dli/li. The sohd hne 
means the process from the strong to the weak couphng phase and the dashed vice versa. 
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Figure 16: Hysteresis for the average hnk-length (/) under the scale-invariant measure Ylidh/h- The 
solid line means the process from the strong to the weak coupling phase and the dashed vice versa. 
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Next, we measured the average link-length (l) as shown in Fig. 1161 where we observe another large 
hysteresis. In the strong-coupling phase, the strong gravitation prevents spacetime manifolds from 
extending widely, resulting in the small value of the length (/). In the weak-coupling phase, conversely, 
spacetime manifolds tend to stretch as widely as possible, resulting in spiky configurations. Such spikes 
of spacetime manifolds often appear also in the weak-coupling phase of quantum RC jl3j . In this sense, 
our data of Fig. is consistent with those obtained in the numerical studies of quantum RC. 

Here we make a short remark on the statistics of our data. To get each point in Figs. El and 1161 
we accumulate 200 configurations, each of which is obtained at an interval of 800 sweeps. Statistical 
errors are also included in the figures, although they are too small to be visible. 

5.2 Fractal structure of 3D pure gravity 

We further investigate other structures of spacetime manifolds under the scale-invariant measure. 
Among them, a fractal structure based on a geodesic distance will be the most important, and its 
validity has been well established in two-dimensional DT [ll] • Here we first explain the fractal structure 
and related issues for our purpose. 

In the case of 3D gravity, the fractal structure is characterized by the surface area distribution 
(SAD) function p{S,D) j45j and the fractal dimension df, their definitions will be given right now 
and closely related to each other. The SAD function is defined as follows |45j : (i) Let us consider a 
connected graph T dual to a three-dimensional closed simplicial manifold T; (ii) The geodesic distance 
between two points in T is defined as the minimum number of steps between the two points; (iii) We 
pick up a point P in the graph T (namely, a 3-simplex in T) and find all points that have a geodesic 
distance D from the starting point P; (iv) The boundary manifold appearing in slicing T at the geodesic 
distance D consists of closed surfaces of various topologies, as shown in Fig. El (v) Then, the SAD 
function p{S, D) is defined as 

p{S, D) = the number of closed surfaces of area S at the geodesic distance D. (38) 

The fractal structure is encoded in the scaling behavior of the SAD function p{S, D) with a scaling 
variable x = S/D", where the exponent a has different values in different phases. p{S,D) x is 
expected to be a function only of x in the scaling regions, and, therefore, it suggests a scaling law 
S ~ D". As a result, the fractal dimension, dj, of spacetime is determined as 

df = a + l . (39) 

Intuitively, such a scaling means a self-similarity of 'time slices' that arise in cutting spacetime manifolds 

at different geodesic distances^^. 

Incidentally, the scaling behavior of p{S, D) was observed numerically in three-dimensional DT !45j, 

while its validity in studying quantum spacetime was first established in two-dimensional DT ,44,. 

Moreover, a similar scaling behavior was reported also in four-dimensional DT j46j . 

Actually, one can think of the geodesic distance D as a time coordinate, corresponding to the so-called temporal 
gauge |47| in 2D quantum gravity. 
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Figure 17: Schematic picture of a three-dimensional fractal manifold that appears in the strong-coupling 
phase of DT. On the 'time slice' at a geodesic distance D one can find many closed surfaces as the 
boundary. The fractal structure is revealed as the scaling property of the distribution of the closed 
surfaces at various geodesic distances. 

According to the values of the gravitational constant G, we can classify the fractal structure into 
three regions: the strong-coupling, the critical, and the weak-coupling regions. An example of the 
fractal manifold with the scaling behavior is schematically depicted on Fig. 1171 which occurs typically 
in the strong-coupling phase of three-dimensional DT. A closed surface of fairly large area appears only 
once at each 'time slice' D, which is the boundary of the so-called mother universe. Concurrently, there 
are many surfaces of small areas at each 'time slice' and they are boundaries of the baby universes. 
The boundaries of the mother and the baby universes consist of closed surfaces of various topologies. 

Having reviewed the fractal structure of lattice gravity, we will devote our attention to the investi- 
gation of the structure in our hybrid model. 

5.2.1 Fractal structure in the strong-coupling limit G = oo 

First, we measure the SAD function p{S, D) in the strong-coupling limit G = cxd in DRC. Numerical 
data obtained is shown in Fig. 1181 where the vertical axis means the SAD scaling function of the form 
p{S,D) X Z)^'^ and the horizontal a scaling variable x = S/D^'^. The different curves show several 
data measured at different geodesic distances D = 5,Q,7 and 8. The boundary surface of the mother 
universe, which has the largest area among all boundaries, has a good scaling property with the scaling 
parameter x = S/D^-^. Accordingly, the fractal dimension is obtained as df = 3.5 + 1 = 4.5. One 
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Figure 18: SAD function p{S, D) in the strong-coupling limit G = oo of 3D pure gravity based on DRC. 
The number of 3-simplices are almost fixed at about = 5000 and each link-length li is integrated 
over a range 1 < li < 10. Different curves correspond to different geodesic distances D = 5,6,7 and 8. 
The scale-invariant measure dk / U is used. 
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can clearly observe in Fig. that the distribution of the mother universe has the Gaussian form. On 
the contrary, the distributions of the boundaries of the baby universes, whose areas are fairly smaller 
than that of the mother, do not exhibit a scaling behavior. In other words, it is impossible that the 
baby part scales in the same way as the mother does; such phenomenon of the existence of two scaling 
variables in a simplicial manifold typically appear also in three-dimensional DT |45j . 

Next, we measure the Euler number of the boundary surfaces. On each closed surface, the Euler 
number x is defined in the usual way: 

X = no - ni + n2 , (40) 

where Uk mean the numbers of fc-simplices {k = 0, 1, 2) on each boundary surface. Our numerical data 
of X sre shown in Fig. 1191 obviously, surfaces with large negative X) exhibiting complicated topologies, 
are identified with the boundary of the mother universe. The distributions of the mother are in the 
shape of mountain (Gaussian) and, in contrast, those of the baby universes have a sharp peak at x = 2. 
This result is consistent with that obtained in the strong-coupling phase of three-dimensional DT j45j . 

How can we draw a physical picture of this region? A plausible answer is a 'confinement of space- 
time', into which all the volume of quantum spacetime is completely confined^'^. Actually, the large 
fractal dimension, whose value is df ~ 3.5 + 1 in this limit, will prevent spacetime from extending, 
although the volume is not small. In addition, the large negative values of x frequently appear in the 
distributions of the mother; this complexity of the topology reminds us of the spacetime foam. 

5.2.2 Fractal structure in the critical region G ~ Gc 

Second, we measure the SAD function and the Euler number distribution at l/lGvrG = 0.6 in the 
critical region. Our numerical data are shown in Fig. 1201 and Fig. 1211 In Fig. I20( the vertical axis is 
the scaling function of the form p{S, D) x D^'^ and the horizontal the scaling variable x = S/D"^'^; the 
exponent 2.3 is vary small compared with the value obtained in the strong-coupling limit. Thus, we 
obtain a smaller fractal dimension df = 2.3 + 1 = 3.3 in this region. In general, the smaller the fractal 
dimension becomes, the more widely the spacetime manifold extends. 

In this region, one cannot separate the contribution of the mother from that of the baby universes, 
because the scaling function p{S, D) distributes smoothly from small x to large one. We have learned 
in two-dimensional DT that the distribution of the mother boundary is universal and, hence, it has 
the close connection with the continuum limit [H]. On the other hand, the distribution of the baby 
boundaries is non-universal; interestingly, the result shown in Fig. [201 is very similar to that obtained 
in two-dimensional DT j44j . although the first-order nature of the transition indicates no continuum 
limit in our case. 

Next, we measure the x distributions of closed surfaces that appear in cutting the simplicial lattice 

at each geodesic distance. Our data is shown in Fig. I2H where the Euler number x is measured at 

several distances. One can see that the x distributions are fairly smoother than those of Fig. 1191 As a 

^^If one can obtain the 'confinement' picture also in four dimensions, it would be regarded as the confinement of 
gravitons. However, we are now in three dimensions, where no gravitons exist. 
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scale-invariant measure Y[idk/k is used. 
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result, we get a milder picture of spacetime in this region and the shape of the x distribution is very 
similar to that obtained in three-dimensional DT j45j . 

5.2.3 Fractal structure in the weak-coupling phase G < Gc 

Lastly, in the weak-coupling phase, we attempted to measure the same quantity p{S, D) in the same 
manner as in the strong-coupling and the critical regions. However, we cannot find any contributions 
from the mother universe in this phase and, what is worse, lattice configurations seem to be almost 
frozen, which have a lower fractal dimension; this phenomenon strongly suggests the spiky nature of 
spacetime in this phase, and such a spike is far from the usual notion of physical spacetime. As is well 
known with numerical studies of DT, this phase corresponds to an elongated branched-polymer with 
a small fractal dimension (for example, df = 2) Similarly, such a spike has been observed also in 
Monte-Carlo studies of quantum RC 13-. 

6 Conclusions and discussion 

We have proposed dynamical Regge calculus (DRC) as a hybrid model of simplicial quantum gravity. 
This model is intended to make physical degrees of freedom larger than those of quantum Regge 
calculus (RC) and dynamical triangulations (DT). Furthermore, the extended model of DRC gives the 
possibility of describing the topology-changing processes of Euclidean spacetime in a dynamical way, 
although there are some difficulties in simulating the processes numerically. 

Algorithmically, the path integral for DRC can be performed through the hybrid (p, q) moves, 
which are an extension of the ergodic (p, q) moves of DT. In particular, the lattice diffeomorphisms 
are generated by the invariance (p, q) moves. It is also an interesting problem that other constructive 
approaches to quantum gravity reproduce the same structure of the lattice diffeomorphisms as that of 
DRC, if one would believe the universality of field theory also in quantum gravity. 

As an application of the lattice diffeomorphism, we tried a lattice-theoretic derivation of the black 
hole entropy. We identified the total number of quantum fluctuations around the event horizon with the 
black hole entropy; such hybrid moves that keep invariant the horizontal geometry play the important 
role. In order to avoid the divergence of the coefficient of the entropy, we required the lattice cutoff /min 
to remain a finite value of order 0(£p). The introduction of the minimal length seems to be consistent 
with a simple interpretation of the space-time uncertainty principle of string theory. 

Moreover, we have carried out the numerical simulations of 3D pure gravity using the two kinds of 
the integration measures. In case of the uniform measure Y\i dh we observed the divergence behavior of 
the lattice size even though we chose the large values of the cosmological constant. This phenomenon 
indicates the exponential unboundedness for the entropy owing to the lattice diffeomorphisms. In this 
case DRC is not well-defined statistical-mechanically. 

In case of the scale-invariant measure Y\^dli/li, however, no pathological behavior occurred in our 
data. Indeed, we calculated the average curvature (R) and the average link-length (l); two pieces of large 
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hysteresis were observed, indicating the existence of the first-order transition between the two phases. 
In the strong-couphng hmit, spacetime manifolds are crumpled to a singular configuration with a large 
fractal dimension. Physically, such a crumpled space might suggest the existence of 'confinement', into 
which spacetime itself is confined. In addition, various topologies appear on time slices of 'confined 
spacetime'. In the critical region, we have obtained the smooth SAD function, resulting in a milder 
picture of spacetime than in the strong-coupling limit. In this case spacetime looks a fractal manifold of 
lower fractal dimension. On the other hand, in the weak coupling phase simplicial configurations become 
spiky, and this phenomenon is essentially the same as that has been observed in the weak-coupling 
phase of quantum RC. Taking into account these theoretical and numerical studies, we conclude that 
DRC can reproduce the numerical results consistent with those of both DT and quantum RC in each 
region. 

Our hybrid model of lattice quantum gravity offers a practical way of studying quantum black 
hole physics and the topology-change of Euclidean spacetime on the lattice. Incidentally, the close 
relation between 3D quantum gravity and quantization of the membrane theory ^ is an interesting 
theme in connection with string theory, though such an attempt generally gives rise to an instability 
problem 59^. Our chief concern is to investigate numerically the physics of strong-coupling quantum 
gravity within the framework of (extended) DRC. 
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